Decoherence in Two-Dimensional Quantum Random Walks with Traps 
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Quantum random walk in a two-dimensional lattice with randomly distributed traps is inves- 
tigated. Distributions of quantum walkers are evaluated dynamically for the cases of Hadamard, 
Fourier, and Grover coins, and quantum to classical transition is examined as a function of the 
density of the traps. It is shown that traps act as a serious and additional source of quantum deco- 
herence. Furthermore, non-trivial temporal dependence of the standard deviation of the probability 
distribution of the walker is found when the trapping imperfections are introduced. 



I. INTRODUCTION 

A simple stochastic process can be introduced by the 
motion of a particle that can move in certain directions 
with some probabilities such that progress of the par- 
ticle is independent of its preceding movements. Find- 
ing the spatial probability distribution of the particle, 
starting such a randomized motion from a given loca- 
tion, is the statement of the so-called random walk prob- 
lem. Random walk and its more generalized extensions 
have been used as simple models for many physical sys- 
tems, in particular in solid-state physics and in astron- 
omy [l[, and in polymer models [2}- Classical random 
walks have found variety of applications in other fields, 
such as in economics and in computational sciences, as 
well. In particular, powerful randomized algorithms, es- 
pecially for graph connectivity, satisfiability, probability 
amplification, [3( and Markovian chain simulations [4j, 
have been developed based on classical random walks. 

Quantum computation and quantum information are 
one of the most active research fields nowadays A 
quantum computer would be capable to run quantum al- 
gorithms, such as Deutsch-Jozsa @, Shor's 0] or Graver's 
algorithm [|[, much faster than any classical algorithms 
running on classical computers. One direction of re- 
search aiming to develop more quantum algorithms is 
to implement some classical algorithms directly by con- 
sidering their quantum analogs. In this respect, quan- 
tum analogs of random walks, the so-called the quantum 
random walks (QRWs) have received much attention re- 
cently (see. e.g., Refs. [1, US EH)- Some new quantum 
algorithms based on QRWs have already been proposed 
0, El, EH EH ■ K is proven that a discrete time quantum 
walk can be used for searching unsorted database with a 
quadratic speedup (l3| . while the continuous time quan- 
tum walks can be exploited for traversing certain graphs 
exponentially faster than any classical algorithm [l2j . 

There are various schemes and systems including ion 
traps [IB] , optical lattices [13, EH , cavity QED [l9l ] , op- 
tical cavity [2(J, and linear optics [U [22j that have 
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been considered for practical realization of QRWs. Using 
nuclear-magnetic resonance two- and three-qubit quan- 
tum information processors, QRW has been demon- 
strated on a square for its continuous time [23[ and dis- 
crete time [13| versions. Such direct implementations of 
QRW are not essential for quantum algorithm develop- 
ments. As long as there exists a quantum computer, 
QRW and associated quantum algorithm can always be 
realized on it. On the other hand, QRWs provide deep 
insight into transition to classical behavior out of a quan- 
tum behavior as well as into decoherence and coherent 
control of quantum systems. Rich quantum dynamics 
available within QRWs makes them appealing systems 
per se. A well-known property of a quantum walker is 
that it is more de-localized than a classical walker such 
that a quantum walker spreads out quadratically faster 
in time than a classical walker. 

Quantum walkers are highly delocalized particles, and 
their QRW is very sensitive to decoherence. The rapid 
transition of QRW to a classical walk may limit their 
implementations and applications, yet at the same time, 
small amount of decoherence can be beneficial to speed 
up quantum algorithms based on QRWs 25]. In order 
to develop optimized quantum walks on noisy quantum 
channels, it is necessary to investigate decoherence mech- 
anisms and characterize their effects on the propagation 
of a quantum walker. Decoherence in QRW can be dis- 
cussed in general terms (for a review, see Ref. [26[) as 
well as in more specialized contexts, depending on a par- 
ticular implementation. For that aim, QRWs on a line 
under various decoherence mechanisms, such as random 
measurements or broken links [27| . and under various 
environmental noises, such as phase flip, bit flip, and 
generalized amplitude damping channels [28] , have been 
examined. The general purpose of these works is to un- 
derstand and resolve the different decoherence problems 
to be encountered in various one-dimensional implemen- 
tations of QRWs. In this work we contribute the same 
purpose by considering an additional decoherence mech- 
anism and by focusing on different QRWs in two dimen- 
sions. Keeping the walkers on the lines or planes of inter- 
est can be nontrivial and practical realizations of QRWs 
may face a serious problem of loss of walkers during the 
evolution. In this work, we consider loss of walkers as a 
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decoherence mechanism. 

A natural imperfection that can cause nonunitary evo- 
lution of a quantum walker would be loss of the walker 
at certain places along the propagation directions. This 
may correspond to absorption of photon in optical im- 
plementations, loss of atom due to thermal fluctuations 
or collisions in atom traps or in optical lattices, or im- 
mobilization of quasi-particles in solids and so on. In- 
deed, in the case of classical RW problem, the possi- 
bility of the trapping of classical random walks by ran- 
domly distributed traps on a lattice has been extensively 
studied using both theoretical and Monte Carlo meth- 
ods [U, M, \M, H [H, [H. The exact source of the 
loss of the quantum walker would depend on a particu- 
lar implementation. In certain systems, such as optical 
lattices, imperfections can be introduced on purpose and 
controllably such as by adjusting the depths of the atom 
traps or by optical addressing of a set of lattice sites. 
The results to be reported in this paper are a systematic 
characterization of effects of immobilizing centers on the 
QRW and are independent of a specific implementation. 

QRW in higher dimensions, especi ally in two dimen- 
sions, have been studied extensively [3a , [37j . In higher 
dimensions, decoherence effect can be of different sig- 
nificance than in one dimension [38], l39j . In this pa- 
per, we numerically investigate the effect of decoher- 
ence that causes quantum to classical transition in a 
two-dimensional (2D) random walk with traps. Typi- 
cal quantum walks, namely, Hadamard (H), Fourier (F) 
and Grover (G) walks are considered. Their relative en- 
durances against decoherence are characterized. Traps 
along the paths of the walkers are recognized as a seri- 
ous source of decoherence that should be reckoned with 
for a potential implementation of the quantum walk, in 
addition to broken links or random measurements. 

The paper is organized as follows. We first present a 
short review of the theory of the two-dimensional QRW 
together with our model of walk on a lattice with traps 
in Sec. HU Results of the numerical simulations and their 
discussions are given in Sec. IIIII Finally, in Sec. IIV1 we 
summarize the conclusions of this work. 



II. QUANTUM RANDOM WALKS IN A 2D 
LATTICE WITH TRAPS 



We consider a QRW in a two-dimensional lattice with 
traps. The traps at a given density are randomly dis- 
tributed before the walk is started then the trap loca- 
tions remain frozen. For a QRW on an infinite two- 
dimensional lattice, the coin space is TL4, spanned by 
the basis states {\j, k) ,j,k S {0,1}}, and the position 
space is 7ioo, spanned by the basis states {\m,n) ,m,n 
integers}. The state of a quantum walker at time t in 
H-4 <S> T~toc space is given by 



One step of the quantum walker is defined by the uni- 
tary operator U such that 



U = S{h ® G), 



(2) 



where I4 is the 4x4 identity matrix, C is the coin oper- 
ator which written as 



c= E E c jkfk >\j,k)\j',k') 

j,k=0 j',fe'=0 



(3) 



and S is the shift operator described by 

S \j, k) \m, n) = \j, k)\m+ {-\y ,n + (-l) fe ) . (4) 

We analyze the results for three different coins which 
are Hadamard, Fourier, and Grover. H, F, and G coins 
in two dimensions are, respectively, given by 



a 



H 



C,: 



C G 




(5) 



(6) 



(7) 



The quantum walker begin to walk at (0, 0) point 
in the two- dimensional lattice. The initial states for 
Hadamard, Fourier, and Grover walks were chosen to 
guarantee a maximum spreading when the walk starts 
at the origin |4(J, and as such for the H, the F, and the 
G coins, they are, respectively, taken to be as 



\m) H 
\m) F 

IV>(o)) G 



~(|00> + i|01)-t|10) + |ll))|0,0), (8) 



((100) 



V2 



|oi) + |io) 



V2 



|11>) |0,0», 



(9) 



-((|00)-|01)-|10) + |11))|0,0)). (10) 



Applying U on the state of the walker ((T|), the rela- 
tion that determines the spatiotemporal evolution of the 
quantum walker is obtained to be as 

1 

Ajkmn (t + 1) = ^ ' Cjkj'k' Ajii;' rnjnj (t), (11) 

j',k'=a 



M*)> = E E A jkmn (t) \j, k) \m,n) . (1) with m 3 = m - (-1)* and n, = n - (-l) fe . With the 



j,k=0 n,m=— c 



amplitude of new positions of the quantum walker given 
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by Eg. ITT|) . the probability distribution of the walker at 
position \m,n) at time t is determined by 

1 

Pm,n{t) = l4/,fc;m,n(*)| 2 - (12) 

3,fc=0 

In the calculations, we assume that the traps are com- 
pletely absorbing, so that if the quantum walker falls into 
a trap, it is annihilated with vanishing probability ampli- 
tudes at any later time. In the next section, the results 
of probability distributions obtained by numerical simu- 
lations for the cases of Hadamard, Fourier, and Grover 
walks with different trap densities will be presented. In 
particular, the examination of the standard deviation a 
will be used to characterize quantum to classical transi- 
tion of the walks. In the case of a classical walk, a oc \/t, 
while a quadratic gain in the spread is achieved in a quan- 
tum walk for which a oc t due to quantum coherence. 
It is known that unitary noises such as broken links or 
nonunitary disturbances such as random measurements 
destroy this gain and make the quantum walk a classical 
one. Here, we explore how the traps, causing nonunitary 



evolution of the quantum walker in the two-dimensional 
lattice, might suppress the benefit of the quantum walk. 
As the traps are randomly distributed, it is necessary to 
consider many different initial configurations of the traps, 
so that we calculate 

1 M 

w = ^E ff - ( 13 ) 

r=l 

where M denotes the number of different configurations 
of random traps. 

III. RESULTS AND DISCUSSIONS 

In our simulations we have used an ensemble of 250 
different configurations of randomly placed traps and the 
ensemble averaged results are presented in the figures for 
the calculated variance. We first investigate the case of 
a Hadamard walk. Quantum walker starts at the ori- 
gin and its walk is numerically simulated for various trap 
densities. Figure [1] shows the probability distribution of 
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FIG. 3: (Color online) Time dependence of the standard devi- 
ation of the Hadamard walk in 2D trapping lattice with trap 
densities p = 0, p = 0.01, p = 0.1, p = 0.25, p = 0.5, and 
classical random walk in 2D lattice. 



the walker at time t = 100 for a trap density of p = 0.01. 
The walk is seen to be still retaining its quantum char- 
acter. One would expect a symmetric centrally peaked 
Gaussian-like distribution about the origin in the classi- 
cal walk. We have found that such a behavior emerges 
and the walk becomes classical for about p = 0.1 at 
t = 100, which is shown in Fig. [51 

The time dependence of the standard deviation for the 
Hadamard walk for different trap densities is plotted in 
Fig. [3] with a log-log scale. When there are no traps, 
p = 0, the highest spatial spread of the walkers or the 
most quantum behavior is available. When trap densi- 
ties increase, the standard deviation is dramatically al- 
tered. For small amount of traps, the spatial spread is 
gradually decreased but still a remains larger than the 
classical one, a c i, and grows faster. However, even when 
a modest amount of traps are introduced, the temporal 
behavior is no longer linear and a can grow as slow as - 
or even slower - than cr c /. This is in stark contrast with 
the other types of decoherence sources, such as broken 
links [1^], for which simple power laws describe the time 
dependence of a. The emergence of the nonlinear be- 
havior of a in a trapped lattice can be attributed to the 
fact that the walkers can survive in a trapped lattice only 
for a time and their probability of survival is related to 
the Kohlrausch- Williams- Watts functions (stretched ex- 
ponentials) of the form exp {-t/rf with < (3 < 1 as 
the stretching exponent and r as a characteristic relax- 
ation time [4j], [42, |43[ . 

When most of the lattice is trapped randomly, so that 
p > 0.5, the walkers become highly localized with a be- 
ing significantly smaller than er c /. In other sources of 
decoherence, when random events are as frequent, a can 
become less than a c i, too [39(. However, traps force a 




FIG. 4: The probability distribution of the Fourier walk for 
p = 0.01 after 100 iterations. 



more severe reduction which increases in time. Although 
a can be smaller in value than the classical one, it cannot 
grow slower than the classical walk in broken link distur- 
bances. With the traps, however, both the magnitude of 
a and its growth rate can be made smaller than those of 
the classical walk. By examining the slopes of the curves 
in Fig. [31 it can be deduced that the transition from 
the quantum to the classical behavior happens at a time 
Tdecoh ~ 5/p, which is the decoherence time. After that, 
the growth rate of a gradually becomes smaller than the 
classical one, which is a reflection of the relatively short 
survival time and a stretched exponential behavior of the 
survival probability. In the broken links caused decoher- 
ence Tdecoh ~ 3/p is found [39[ . This suggests that two- 
dimensional lattices can endure only slightly longer to 
decoherence due to traps. 

Let us now look at the other typical quantum walks. 
The probability distribution of the walker in the Fourier 
walk when the trap density is p = 0.01 is given in Fig. 0] 
at iteration t = 100. Analogous to the Hadamard walk, 
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FIG. 5: The probability distribution of the Fourier walk for 
p = 0.1 after 100 iterations. 



a significant quantum behavior is maintained at such low 
trap densities, while as shown in Fig. [5] for p = 0.1 and 
t = 100 the classical behavior is emerged. The time 
evolution of the standard deviation for the Fourier walk 
for different trap densities is depicted in Fig. [6] with 
a log-log scale. A strong nonlinear behavior is found 
in this case as in the Hadamard walk. When Fig. [5] 
is compared to Fig. [3J it is seen that the Fourier walk 
becomes classical more quickly than the Hadamard walk. 
This observation is in agreement with the broken links 
case [39l | . 

Finally, we look at the case of the Grover walk for 
which Fig. [7] describes the probability distribution at 
time t = 100 for the case of a trap density p = 0.01, and 
Fig. [8] for p = 0.1. The low and the high trap density 
behaviors lead to the similar conclusions as in the Grover 
and the Fourier walks, while the time dependence of the 
standard deviation in Fig. ^indicates that the Hadamard 
walk comes to the classical behavior more slowly than 
the Fourier and the Grover walks. The Hadamard walk 
in two-dimensional lattice can be said to be more durable 
in dccohcrcncc due to traps while the Fourier is the least 
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FIG. 6: (Color online) Time dependence of the standard de- 
viation of the Fourier walk in 2D trapping lattice with trap 
densities p — 0, p = 0.01, p = 0.1, p = 0.25, p = 0.5 and 
classical random walk in 2D lattice. 
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FIG. 7: The probability distribution of the Grover walk for 
p = 0.01 after 100 iterations. 
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FIG. 8: The probability distribution of the Grover walk for 
p = 0.1 after 100 iterations. 



durable. This observation is also valid for the case of 
decoherence due to broken links [39], which is in fact a 
fundamentally based on symmetry properties of the prob- 
ability distributions [39j when there is no decoherence 
and thus should be valid for any type of imperfections 
introduced. 
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FIG. 9: (Color online) Time dependence of the standard de- 
viation of the Grover walk in 2D trapping lattice with trap 
densities p — 0, p = 0.01, p = 0.1, p = 0.25, p = 0.5 and 
classical random walk in 2D lattice. 



the traps as a source of decoherence reveals that trap- 
ping forces transition to a classical walk from a quantum 
walk in a relatively longer time in comparison to the bro- 
ken links. While all three types of the quantum walks 
become classical after a time decreasing with the trap 
density, it is found that quantum walk with a Hadamard 
coin is found to be the most tolerant to the density of 
the traps. A quantum walker with a Fourier coin would 
exhibit the lowest endurance against traps. Major dif- 
ference in a trapped lattice is the survival time of the 
walkers that causes slower growth rates of the standard 
deviations than classical walks and lack of simple power- 
law temporal growth of the standard deviations due to 
the stretched exponential behavior of the survival prob- 
abilities of the walkers. In addition to broken links and 
random measurements, traps corresponding to defects 
and impurities thus should be considered as a serious 
and additional source of decoherence in practical imple- 
mentations of quantum random walks. 



IV. CONCLUSIONS 

Introducing traps along the path of a random walker 
has been investigated in detail classically in the context 
of the survival probability of particles in solid-state sys- 
tems with defects and impurities. We examine this prob- 
lem for the case of quantum random walk and show that 
traps act as a significant source of quantum decoherence 
and make the walk a classical one even when a mod- 
est number of traps are present. We consider the cases 
of three different quantum coins, which are Hadamard, 
Fourier, and Grover. The corresponding initial states for 
each coin are chosen to guarantee a maximum spreading 
when the walk starts at the origin. Characterization of 
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